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We propose an importance-truncation scheme for the large-scale nuclear shell model that extends its range 
of applicability to larger valence spaces and mid-shell nuclei. It is based on a perturbative measure for the im¬ 
portance of individual basis states that acts as an additional truncation for the many-body model space in which 
the eigenvalue problem of the Hamiltonian is solved numerically. Through a posteriori extrapolations of all ob¬ 
servables to vanishing importance threshold, the full shell-model results can be recovered. In addition to simple 
threshold extrapolations, we explore extrapolations based on the energy variance. We apply the importance- 
truncated shell model for the study of ^®Ni in the pf valence space and of ®®Zn and ^''Ge in the pfgy, space. We 
demonstrate the efficiency and accuracy of the approach, which pave the way for future shell-model calculations 
in larger valence spaces with valence-space interactions derived in ab initio approaches. 

PACS numbers: 21.60.Cs, 21.10.-k, 27.50.-He 


Introduction. The nuclear valence-space shell model is 
one of the work horses in nuclear structure theory. It is very 
successful for the description of spectra and spectroscopic ob¬ 
servables over a large range of nuclei and plays an important 
role in guiding and interpreting experiments from stable to 
exotic nuclei lHHIll. Two critical aspects in the application 
of the shell model (SM) are the construction of the effective 
valence-space interaction as well as corresponding effective 
operators and the solution of the eigenvalue problem in the 
model space of the valence nucleons. 

Traditionally, effective valence-space interactions are con¬ 
structed using renormalized nucleon-nucleon interactions 
combined with phenomenological fits of matrix elements to 
nuclei within the valence space lH] [3l [HI [HI . Though this 
phenomenological approach allows for a rather accurate de¬ 
scription, it lacks a rigorous connection to the underlying 
nuclear interaction and does not provide a consistent frame¬ 
work for the treatment of observables other than the energy. 
Recently, a set of novel approaches to systematically derive 
valence-space interactions and operators have emerged m- 
l23l . They offer new insights into valence-space interactions 
and can be linked to ab initio calculations. 

Once the valence-space interaction is specified, the SM re¬ 
duces to the solution of a large-scale matrix eigenvalue prob¬ 
lem. Its dimension grows combinatorially with the number of 
valence orbitals and nucleons. Starting with valence spaces 
covering the p/-shell, the m-scheme model spaces reach di¬ 
mensions beyond 10® around mid-shell, which is approach¬ 
ing the limits of present computational approaches for sparse 
eigenvalue problems. When going to valence spaces covering 
more than one major shell, the model-space dimension poses a 
severe limitation to the applicability of the SM. In these cases 
additional truncations or more sophisticated methods like the 
Monte Carlo shell model (MCSM) Il24l 1^ or a density ma¬ 
trix renormalization group treatment of the SM li2^ have to 
be employed. 

We propose the importance-truncated shell model (IT-SM) 
to overcome this limitation. It combines the SM with an 
importance-truncation scheme that is successfully applied in 


no-core configuration-interaction approaches for some time 
ETII . In addition we use refined extrapolation schemes based 
on the energy variance to reduce the uncertainties of the IT- 
SM calculations. Together, importance truncation and ex¬ 
trapolation provide an accurate tool for systems and valence 
spaces beyond the reach of standard SM calculations. 

Importance Truncation. The importance truncation is a 
physics-driven, adaptive truncation of the many-body model 
space based on a measure for the importance of individual 
basis states for the description of a specific set of eigen¬ 
states of a given Hamiltonian. The importance measure is de¬ 
fined through the amplitude of the individual basis states in 
the expansion of the eigenstates, obtained a priori in lowest- 
order many-body perturbation theory. By imposing a thresh¬ 
old with respect to this importance measure we define an 
importance-truncated model space tailored specifically for 
the target eigenstates and Hamiltonian under consideration. 
Eventually, variations of the importance threshold and extrap¬ 
olations to vanishing threshold can be used to extract observ¬ 
ables in the limit of the full model space. This scheme is 
applied very successfully in the context of the no-core shell 
model (NCSM) Il27ll28l . 

The construction of the importance-truncated space is 
based on a set of reference states which are obtained 

from a previous diagonalization in a small space, that repre¬ 
sent the target eigenstates. The basis states that contribute 
to the reference states |'I'^™f) span the reference space AJref- 
We estimate the importance of basis states |®y) outside Ai^ef 
by means of the amplitudes Ky'^ - - {®y| H |'I'|”j )/Acy of 
the first-order perturbative correction to I'Tref), where the en¬ 
ergy denominator Aey corresponds to the unperturbed single- 
particle excitation energy of the basis state |®y). Only basis 
states with importance measure |/Cy"^| larger than a given im¬ 
portance threshold /c^in for at least one reference state |'I'|™^) 
are included in the importance-truncated model space. 

In the case of a two-body Hamiltonian, the simple first- 
order importance measure cannot probe basis states that differ 
by more than a two-particle-two-hole (2p2h) excitation from 
any state in M^ef- Therefore, we embed the construction of 
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the importance-truncated space into an iterative scheme. For 
the valence-space SM, we use the number of valence parti¬ 
cles above the orbits that are (partially) occupied in the lowest 
energy configurations to define a truncation parameter r^ax- 
For Fmax = 0, all Slater determinants with valence nucleons 
distributed in the lowest accessible orbits are in the model 
space. For Fmax = 2, up to two valence nucleons are pro¬ 
moted to higher-lying orbits—this model space can be gener¬ 
ated through Iplh and 2p2h excitations on top of the T^ax = 0 
space. Thus, we combine a sequential increase of the tmnca- 
tion parameter r^ax with the importance-selection procedure. 
This sequential IT-SM scheme is analogous to the sequential 
IT-NCSM scheme discussed in Ref. lIZTl . 

The complete IT-SM calculation proceeds as follows: We 
start with a conventional SM calculation for small Tmax, e.g. 
Tmax = 0, and select a set of target eigenstates. We define the 
reference states ) by filtering the important components 
of these eigenstates through a so-called reference threshold 
Cmin with respect to the amplitudes from the SM calculation. 
With these reference states we construct importance-truncated 
spaces with Tmax = 2 for a sequence of importance thresholds 
ATmin- In each space we solve the eigenvalue problem and com¬ 
pute the relevant observables. The eigenvectors for the largest 
importance-trancated space define the new reference states, 
again imposing a reference threshold Cmin, for constructing 
the importance-trancated spaces for - 4. This proce¬ 
dure can be iterated until Tmax reaches the number of valence 
particles and thus probes the full model space. In the limit 
(/Cmin, Cmin) —> 0, tWs algorithm is guaranteed to reproduce 
the results in the full model space at each Tmax- 

The results of IT-SM calculations for different thresholds 
/Cmin and Cmin are depicted in Fig.[2 As a test case, we consider 
^®Ni in a pf valence space using the gxpfIa interaction 123 . 
The full /«-scheme dimension of this model space is 1.09x10®, 
which is at the limit of routine SM calculations. The results 
presented in Fig. [T] show the dimensions and the lowest en¬ 
ergy eigenvalues as function of ATmin- Note that the energy axis 
is extremely magnified and spans only 80 keV. The dimen¬ 
sions of the importance-truncated spaces are reduced drasti¬ 
cally, by about two orders of magnitude as compared to the 
full SM space. At the same time, the absolute energies in 
the largest importance-truncated spaces, corresponding to the 
smallest /Cmin and Cmin thresholds, differ by only about 10 keV 
from the full SM. This demonstrates the efficiency of the im¬ 
portance truncation—it separates the 10® basis states that de¬ 
termine the bulk of the energy from the 10® basis states that 
are responsible for the residual 10 keV. 

Threshold Extrapolation. We can approximately account 
for the effects of basis configurations excluded from the 
importance-trancated spaces by an a posteriori extrapolation 
of the observables. The simplest extrapolation addresses the 
importance threshold /tmin- Since the energy eigenvalues de¬ 
pend smoothly on /fmin we can fit simple functions to the set of 
energies obtained for different /Cmin values and extract the en¬ 
ergies for /fmin —> 0. Since we do not have a theoretical model 
for the functional dependence on the importance thresholds. 
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FIG. 1. (color online) Dimension of the importance-truncated model 
space (a) and ground-state energy relative to the core (b) for ^®Ni in 
the pf valence space as a function of the importance threshold for 

reference thresholds Cmin = {!(•), 2 ( 4 ), 3 (^fjx 10“® and Tmax = 16 

using the gxpfIa interaction. The model space has been constructed 
for the simultaneous description of the six lowest eigenstates. For the 
threshold extrapolation we use polynomials of order two and three. 
The red lines denote the full ///-scheme dimension and the ground- 
state energy of the full SM (2). 


we use simple polynomials, typically of order two to four. 
Varying the order of the polynomials gives an estimate for 
the uncertainty of this threshold extrapolation. In Fig.[TJb) we 
have included examples for fits with second and third-order 
polynomials for the ground-state energies of ^®Ni. Note that 
the uncertainty of the K^in extrapolation is small compared to 
the residual dependence on the reference threshold Cmin- 

This simple threshold extrapolation does not require addi¬ 
tional computations and can be applied to all observables on 
equal footing (cf. Fig. [^. However, it exclusively addresses 
the importance threshold /Cmin and uncertainties of the poly¬ 
nomial extrapolations can be sizeable. One can improve on 
this by including additional information on the excluded ba¬ 
sis states, e.g., through a second-order perturbative estimate 
of their contribution to the energy, as done successfully in the 
IT-NCSM (see Ref. Ell for details). 

Variance Extrapolation. As an alternative to the simple 
threshold extrapolation, we consider a more elaborate ex¬ 
trapolation scheme based on the energy variance AE^ = 
^vj/|h2|vj/) _ (vf'iHI'T)®, which was used in the SM context be¬ 
fore II 30 I 432 I . By construction, the energy variance vanishes 
for the exact eigenstates and, thus, serves as a measure for 
the distance of an approximate state obtained in a truncated 
subspace from the energy eigenstate in the full space. As dis¬ 
cussed in Ref. 021, the energy is expected to show a predom¬ 
inantly linear dependence on the energy variance, with sub¬ 
leading quadratic corrections. We thus have a simple model 
and a robust two- or three-parameter fit function at hand that 
provides accurate extrapolations. 

The energy variance captures nontrivial information on the 
full model space through the expectation value 
This is seen by inserting an identity operator represented 
in the full model space in between the product of the two 
Hamiltonians—the variance explicitly probes the coupling to 
states outside of the truncated subspace. In practical calcu¬ 
lations we can choose the full space we wish to extrapolate 
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FIG. 2. (color online) Energy-variance extrapolation of the ground- 
state energy relative to the core for ^®Ni obtained in IT-SM using 
the gxpfIa interaction. In panel (a) results for different reference 
thresholds Cmm = {!(•), 2(4), 3 (b)) x lO^"* for F^ax = 16 are 
shown. In panel (b) calculations for different truncations = 
j4(«), 6(4), 8(b), 10(4)1 with Cmin = 1 x 10^“* are depicted. 


to. In the context of the IT-SM it is desirable to extrapolate 
to the full SM space defined by the valence orbits without ad¬ 
ditional truncations. In this way the variance extrapolation 
remedies all truncations used in the IT-SM, i.e., the A-min, Cmin, 
and Tmax truncations. Therefore, the variance extrapolation is 
much more powerful than the simple threshold extrapolation. 

In Fig. S the variance extrapolation of the ground-state en¬ 
ergy of ^^i is illustrated, where panel (a) shows se¬ 
quences for different reference thresholds Cmin and panel (b) 
shows Amin Sequences for different Tmax truncations. The hrst 
remarkable observation is that the Amin sequences for different 
Cmin fall onto a straight line. Consequently the variance ex¬ 
trapolations for the different Cmin give the same result. The 
variance-extrapolated energy is in excellent agreement with 
the result for the full space reported in Ref. . Even with an 
additional Tmax truncation, as shown in Fig. |^b), the results 
beyond Tmax = 6 fall onto the same line. For severe trunca¬ 
tions, e.g. Tmax = 4, we observe larger energy variances that 
cannot be extrapolated reliably. 

The advantages of the variance extrapolation are that a sim¬ 
ple and robust ht model is available and that the extrapolation 
remedies all truncations inherent to an IT-SM calculation. The 
disadvantage is that substantial computational effort goes into 
the evaluation of the energy variance, typically the compu¬ 
tation of the variance needs more computing time than the 
complete IT-SM calculation. 

Applications. Using the IT-SM with threshold and vari¬ 
ance extrapolation we now discuss the spectroscopy of ^®Ni in 
the pf shell with the gxpfIa interaction lisa. We demonstrate 
the robustness of the IT-SM by comparing energies and elec¬ 
tromagnetic observables with full SM results obtained with 
the Antoine code |[I]|34l|35l or extracted from Ref. ||2l . 

Figure shows the excitation spectrum of ^®Ni for three 
different reference thresholds as function of energy variance 
and importance threshold. On the scale of typical excitation 
energies the Amin and Cmin dependence is very weak. Both, 
the variance- and the threshold-extrapolated energies are in 



FIG. 3. (color online) Energy-variance (a) and threshold (b) extrap¬ 
olation of the energies of the six lowest natural-parity states of ^^Ni 
using the gxpfIa interaction with Cmin = {!(•), 2(4), 3(b))x 10-* 
and Tmax = 16. For the variance and threshold extrapolations, poly¬ 
nomials of order two and three have been employed, respectively. 
The red lines show the full SM results extracted from III 



FIG. 4. (color online) Natural-parity spectrum of ^*Ni as a function 
of Tmax in terms of absolute energies relative to the core computed in 
the IT-SM with Cmin = 1 x 10““* using the gxpfIa interaction. The 
right-hand columns show the results of an energy-variance extrapo¬ 
lation (A£^) and the full SM energies extracted from lE). 

excellent agreement with the full SM results where available. 
The second 0^ state, however, shows a quite distinct behavior. 
Its Amin and Cmin dependences are stronger than for all other 
states and the energy variances are signihcantly larger. This 
indicates a particularly complicated structure, in this case due 
to deformation, resulting in many small components in the 
basis expansion of the eigenstate and, thus, a less accurate 
approximation in the importance-truncated space. The sim¬ 
ple threshold extrapolation does not capture the contribution 
of all these small components and cannot correct for the size¬ 
able Cmin dependence. The variance extrapolation, however, 
provides a reliable extrapolation and even restores the correct 
level ordering in excellent agreement with the full SM. Partic¬ 
ularly for these fragile states, the variance extrapolation offers 
signihcant advantages. 
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FIG. 5. (color online) Threshold dependence and extrapolation for 
the quadrupole moment of the 2| state (a) and the S(E2 : 2| —> 
0|) transition strength (b) for ==Ni. The wave functions have been 
obtained in an IT-SM calculation using the gxpfIa interaction for 
Fmax = 8 and the reference thresholds Cn,i„ = {1 (•), 2 (♦), 3 (■)) x 
IQ-**. The red lines represent the full SM results obtained with the 
Antoine code ElMlEll. 


Figurej^summarizes the extrapolated energies for the low¬ 
est six natural parity states of ^®Ni. The results of threshold 
extrapolations for a sequence of Tmax-truncated calculations 
are shown in the main part of the plot, followed by the spec¬ 
trum obtained from the variance extrapolation with Tmax = 8 
and the full SM result 121. Starting from Tmax = 8 the spec¬ 
trum is rather stable and in good agreement with the full SM 
results, except for the second 0=^ state discussed above. The 
energy-variance extrapolation for r^ax = 8 yields excellent 
agreement with the full SM for all states. 

Since the IT-SM also provides the eigenstates in the 
importance-truncated space, we have access to all other ob¬ 
servables, particularly to electromagnetic moments and transi¬ 
tions relevant for spectroscopy. For each we compute the 
observable of interest using the respective eigenvector. Figure 
[^illustrates the dependence of the quadrupole moment and the 
B(E2) transition strength from the first 2^ state to the ground 
state in ^®Ni on the importance threshold. Also these observ¬ 
ables show a smooth dependence on iCmin and allow for simple 
polynomial extrapolations to vanishing importance threshold. 
There is a mild dependence of the A-min-extrapolated results on 
Cmin which is of the same magnitude as the uncertainty of the 
Amin extrapolation. Within these small uncertainties, the ex¬ 
trapolated quadrupole moment and B{E2) transition strength 
are in excellent agreement with full SM calculations proving 
that spectroscopic observables are also directly accessible in 
the IT-SM. We note that extrapolations using the energy vari¬ 
ance do not improve these results. 

We conclude this discussion with a first application of the 
IT-SM in a valence space covering more than one major shell. 
This will be an important future field of application of the 
IT-SM in conjunction with the new valence-space interactions 
derived in ab initio approaches. We consider a p/gy, valence 
space using the pfg9b3 interaction 1361 |37l and study ®°Zn 
and ^'^Ge with full model-space dimensions of 2.2 x lO'^ and 
1.7 X 10'"^, respectively. Since these extended model spaces 
are susceptible to center-of-mass spuriosities, we use a Law¬ 



FIG. 6. (color online) Lowest natural-parity states of ®°Zn (a) and 
=^Ge (b) computed in the IT-SM for Cmin = 2 x IQ-"^ using the pfg9b3 
interaction with subsequent threshold extrapolation for different val¬ 
ues of Tmax- The right-hand columns show the results obtained from 
the energy-variance extrapolation (AL^) of the Tmax = 6 results. 
The dashed line shows an approximation for the energy of the 2* 
state calculated from the excitation energy obtained in the IT-SM for 
Tmax = 10 and the -extrapolated ground-state energy. For ^''Ge, 
the MCSM results I31II32I are shown for comparison. 


son prescription to diagnose center-of-mass contaminations 
|l38l. Particularly ®^Ge has been studied before in the MCSM 
Emm using the same interaction. Very recently, a study us¬ 
ing the density-matrix renormalization group with the SM tar¬ 
geted the same nucleus and valence space ll26l . These compet¬ 
ing approaches highlight the difficulty of these calculations. 

Figure [^ shows the absolute energies of the lowest states 
in ®°Zn and ®^Ge extracted from a simple threshold extrap¬ 
olation for a sequence of Tmax truncated spaces and from an 
energy-variance extrapolation. Whereas the spectra seem con¬ 
verged at Tmax = 8 for ®°Zn there is still some dependence on 
Tmax for ^'^Ge. Moreover, for ®"^Ge the variance extrapolation 
gives the ground-state energy about 0.5 MeV lower than the 
threshold-extrapolated energy at Tmax =10, due to effects of 
the Cjnax and Tmax truncations ignored in the threshofd extrap- 
ofations. The sensitivity to these truncations resufts from the 
strong deformation of ®'^Ge, which requires many smaff com¬ 
ponents in the SM basis expansion to describe. The variance 
extrapofation captures these subtie effects and yieids exceiient 
agreement with the MCSM resuits EUEll- 

Conclusions. We have introduced the IT-SM approach 
and demonstrated its ability to extend the reach of valence- 
space SM calculations into the domain of large valence spaces 
and mid-shell nuclei. In addition to the threshold extrapola¬ 
tion, we adopted an extrapolation in terms of the energy vari¬ 
ance for the first time in the IT context. Generally, the thresh¬ 
old extrapolation provides sufficiently accurate energies and 
electromagnetic observables at no extra computational cost. 
In specific cases, e.g. for states governed by deformation, the 
energy-variance extrapolation provides better accuracy for en¬ 
ergies at significant extra cost. 

The IT-SM framework is ideally suited to study valence 
spaces spanning two or more major shells with effective in¬ 
teractions derived in an ab initio framework, such as the in¬ 
medium similarity renomalization group EOl |2TI| or the Lee- 
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Suzuki approach 122] 123]. Together, these new developments 
offer unique perspectives for detailed nuclear structure inves¬ 
tigations beyond the reach of the conventional SM. 
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